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Abstract 

Let ^ n be a product of n independent, identically distributed random matrices 
M, with the properties that ^f n is bounded in n, and that M has a deterministic 
(constant) invariant vector. Assume that the probability of M having only the 
simple eigenvalue 1 on the unit circle does not vanish. We show that $„ is the sum 
of a fluctuating and a decaying process. The latter converges to zero almost surely, 
exponentially fast as n — > oo. The fluctuating part converges in Cesaro mean to a 
limit that is characterized explicitly by the deterministic invariant vector and the 
spectral data of E[M] associated to 1. No additional assumptions are made on the 
matrices M; they may have complex entries and not be invertible. 

We apply our general results to two classes of dynamical systems: inhomoge- 
neous Markov chains with random transition matrices (stochastic matrices), and 
random repeated interaction quantum systems. In both cases, we prove ergodic 
theorems for the dynamics, and we obtain the the limit states. 

1 Introduction 

In this paper we study products of infinitely many independent, identically distributed 
random matrices. The matrices we consider satisfy two basic properties, reflecting 
the fact that they describe the dynamics of random quantum or classical dynamical 
systems. The first property is that the norm of any product of such matrices is bounded 
uniformly in the number of factors. It reflects the fact that the underlying dynamics is 
in a certain sense norm-preserving. The second property is that there is a deterministic 
invariant vector. This represents a normalization of the dynamics. Two important 
examples of systems falling into this category are inhomogeneous Markov chains with 
random transition matrices, i.e. products of random stochastic matrices, as well as 
repeated interaction open quantum systems. In this paper we present general results 
on infinite products of random matrices, and applications of these results to the two 
classes of dynamical systems mentioned. Our main results are convergence theorems 
of the infinite random matrix product, Theorems 1.1, 1.2 and 1.3. They translate into 
ergodic theorems for the corresponding dynamical systems, Theorems 1.4 and 1.5. 
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1.1 General Results 

Let M(co) be a random matrix on C d , with probability space (Q, p). We say that 
M(uj) is a random reduced dynamics operator (RRDO) if 

(1) There exists a norm ||| • ||| on C d such that, for all to, M(uj) is a contraction on 
C d endowed with the norm 1 1 1 • 1 1 . 

(2) There is a vector ijjs, constant in u>, such that M(u>)ips = i/js, for all uo. 

We shall normalize %ps such that ||Y>«s|| = 1 where || • || denotes the euclidean norm. 
To an RRDO M(w), we associate the (iid) random reduced dynamics process (RRDP) 

tf n (w) :=M(wi)---M(w n ), Zoen N \ 

We show that ^> n has a decomposition into an exponentially decaying part and a 
fluctuating part. To identify these parts, we proceed as follows. It follows from (1) 
and (2) that the spectrum of an RRDO M(u) must lie inside the closed complex unit 
disk, and that 1 is an eigenvalue (with eigenvector ips)- Let Pi(w) denote the spectral 
projection of M(uj) corresponding to the eigenvalue 1 (dim Pi (a;) > 1), and let Pi(u) 
be its adjoint operator. Define 

y>H := PiHVs, (l.i) 

and set 

p{u) = IVsWHi. 

For ip, <fi S C d , we denote by the rank-one operator C d B \ \—> \ip)((f>\x = (fii X) V'j 

and our convention is to take the inner products linear in the second factor. We put 

Q(w) = 1-P(w). 

Note that the vector ip(uj) is normalized as (ips,ip(uj)} = 1. We decompose M{u) as 

M(w) = P(w) + Q(w)M(a;)Q(w) =: P(w) + M q (lu). (1.2) 

Taking into account this decomposition, we obtain (c.f. Proposition 2.1) 

* n (aJ) := M(a;i) • • • M{uj n ) = \^ s ){6 n {uj)\ + M Q (a; 1 ) • • • M Q (u n ), (1.3) 

where ^ n (aJ) is the Markov process 

S n (ai)=M*K)-M*( W2 )^ 1 ), (1.4) 

M*(wj) being the adjoint operator of M(ujj). We analyze the two parts in the r.h.s. of 
(1.3) separately. 

Let M.(e) be the set of RRDOs M whose spectrum on the complex unit circle 
consists only of a simple eigenvalue {1}. 

On N * we define the probability measure dP in a standard fashion by 

dP = IL;>idpj, where dpj = dp, Mj G N*. 
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Theorem 1.1 (Decaying process) Let M{uj) be a random reduced dynamics opera- 
tor. Suppose that p(M(w) G M-(e)) > 0. Then there exist a set f2i C S7 N * and constants 
C,a > s.t. P(f2i) = 1 and s.t. for any ZJ G fii, i/iere exists an n (w) s.i. /or any 
n > no(uJ), 

\\M Q (u Jl )---M Q (u; n )\\<Ce- an . (1.5) 

Remarks. 1. In the case where M(u>) = M is constant, and M G M.<e), one readily 
shows that for any e > there is a C t such that ||(M Q ) n || < C e e~ n ^- e \ for all n > 0, 
and where 7 = min-jggpg^^yix} | log \ z\ | (see e.g. [7], Proposition 2.2). It is remarkable 
that in the random case, the mere condition of M having an arbitrarily small, non- 
vanishing probability to be in M.^e) suffices to guarantee the exponential decay of the 
product in (1.5). 

2. Any stochastic matrix whose entries are all nonzero belongs to M^e)- 

3. If {1} is a simple eigenvalue of M{oj) then the decomposition (1.2) is just the 
spectral decomposition of the matrix M{u). 

4. The choice (1.1) ensures that 4>(oj) is an eigenvector of M*{oj). Other choices 
of measurable ip(uj) which are bounded in oj lead to different decompositions of M(w), 
and can be useful as well. For instance, if M(u) is a bistochastic matrix, then one can 
take for ip{oo) an M*(w)-invariant vector which is constant in u. 

Our next result concerns the asymptotics of the Markov process (1.4). Set E[/] = 
f n f(io)dp(u>) for a random variable /, and denote by P\ t E[M] t ne spectral projection 
of E[M] onto the eigenvalue {1}. 

Theorem 1.2 (Fluctuating process) Let M(uj) be a random reduced dynamics op- 
erator. Suppose that p(M((j) G M.(e)) > 0- Then we have E[M] G M(e)- Moreover, 
there exists a set CI2 C N s.t. P($^2) = 1 and, for all uJ G 

1 N 

hm ^0 n( aJ) = 0, (1.6) 

N^oo iv * — * 
n=l 

where 

e = (i- eimq]*)- 1 e[v»] = Pt m] m = Pi nM] i>s. (1.7) 

Remarks. 5. In the case where M is constant in u, we have E[Mq]* = (Mq)*, K[tp] = ip, 
and under the assumption of Theorem 1.2, that M G Mie)- Therefore, P± = P = 
{ips) (V 7 ! an d hence Q*ip = 0, and (Mq)*iP = 0. Consequently, we have 6 = tp. This 
coincides with the results of [7]. As we will see in Section 3.4, the latter equality is not 
satisfied for general, w-dependent matrices M. 

6. The ergodic average limit of n (jD) does not depend on the particular choice of 
ip{uj)- This follows from the last equality in (1.7). 

7. We show in Proposition 2.5 that for every fixed U, 9 n (uj) converges if and only if 
ip(u} n ) converges, and that the limits coincide if they exist. 

Combining Theorems 1.1 and 1.2 we obtain the following result. 
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Theorem 1.3 (Ergodic theorem for RRDP) Let M(lo) be a random reduced dy- 
namics operator. Suppose p(M(w) G M.(e)) > 0. Then there exists a set Q3 C Q w s.t. 
P(0 3 ) = 1 and, for all ZJ £ Q 3 , 

1 N 

\im^-J2 M ("i)--- M (. UJ n) = \ips)(e\ = Pi,m[M]- (1-8) 

n=l 

Remark. 8. If one can choose ip(u) = V to be independent of u> (see also remark 
4 above), then we show (see (2.8)) that 9 n (ZJ) = ip, for all n,ZJ. It thus follows from 
(1.3)-(1.5) that lim^oo M(uj\) ■ ■ ■ M{ui n ) = |V , s)(V'l> a - s - ; exponentially fast. 

1.2 Inhomogeneous random Markov chains 

Products of random matrices satisfying conditions (1) and (2) are important in the 
study of inhomogenous random Markov chains. A finite Markov chain is a system 
consisting of d < 00 states on which a stochastic process is defined in the following 
way. At time step n, n £ N*, the probability to make a transition from state i to state 
j, where i,j £ {1, . . . , d}, is given by pff ■ Hence the relation Y^j=iPij^ = 1- Arranged 

into a d x d matrix (M n ) y = , these probabilities yield a stochastic matrix. The 
probability of jumping from state i to state j between times and n is given by the 
matrix element ij of the product 

V n = M 1 M 2 ---M n . (1.9) 

The Markov chain (1.9) is called inhomogeneous if the transition probabilities depend 
on the time step n. We have that for all n, 

- M n is a contraction for the norm |||M||| = maxi=i t ... j dYlj=i \(^)ij\-> anc ^ 

- the vector ^5 = ^(1, 1, ■ ■ ■ , 1) T is a normalized eigenvector of M n , associated to 
the eigenvalue 1. 

Consequently, if the transition probabilities are chosen randomly at each time step, 
the corresponding operator M{uj) is an RRDO, and the product (1.9) is an RRDP. 
If > 0, for all then the Perron- Frobenius Theorem asserts that 1 is a simple 
eigenvalue of M n , and that all other eigenvalues lie inside the open complex unit disk. 
Thus the conclusions of Theorems 1.1 -1.3 hold for M{oj) a random stochastic matrix, 
which, with positive probability, has strictly positive entries. 

Results on convergence of products of random stochastic matrices, or non-negative 
matrices, are numerous, see e.g. the references in Section 1.4. However, these results 
mostly concern properties of the limiting distribution, if it exists, in terms of the 
properties of the distribution of the RRDOs. Those results rely heavily on the positivity 
of the elements of the considered matrices. When studying convergence in distribution, 
the order of the factors in the product does not matter, and usually products of the 
form 

$ n = M n M„_ 1 ---M 1 (1.10) 
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as studied as well (compare with (1.9)). While products of the form (1.10) are not our 
main concern (and are easier to study), our techniques still yield results for them that 
are stronger than those for products (1.9). We show in Section 3.5 a general result on 
products (1.10) which, applied to products of stochastic matrices, yields the following 
result. 

Theorem 1.4 Suppose that M{u) takes values in the set of stochastic matrices, and 
that p(M(uj) G Ve) > 0, where Ve denotes the set of stochastic matrices with stricly 
positive entries. Then there exists an a > and a set ^4 C S1 N *, such that ¥(0,4) = I, 
and, for all To 6 ^4, 

<D n (uJ) = MK)MK_!) • • • M(wi) = IVttMooMI + 0^(e- Qn ), (1.11) 

where r?oo( w ) a vector valued random variable, where tps = T^O-i 1> " " " > 1) T ; anc ^ 
where the remainder term depends on uJ £ Q4. 

Remark. 9. As expected, Efryoo] = 0, with 6 given in (1.7). 

1.3 Random repeated interaction open quantum systems 

In this section we present the mathematical framework of open quantum systems. We 
establish the link between the dynamics and infinite products of random matrices, and 
we present an application of Theorems 1.1-1.3. 

1.3.1 Definition of repeated interaction quantum systems 

A repeated interaction quantum system consists of a subsystem S which interacts 
successively with the elements £ m of a chain C = £ 1 + £ 2 + ■ ■ ■ of independent quantum 
systems. At each moment in time, S interacts precisely with one £ rn (m increases 
as time does), while the other elements in the chain evolve freely according to their 
intrinsic (uncoupled) dynamics. The complete evolution is described by the intrinsic 
dynamics of S and of all the £ m , plus an interaction between S and £ m , for each m. The 
latter consists of an interaction time r m > 0, and an interaction operator V rn (acting 

on S and £ m ); during the time interval [t± + ■ ■ ■ + T m _i, t\ H + r m ), S is coupled to 

£ m via V m . Systems with this structure are important not only from a mathematical, 
but also from a physical point of view. They arise in fundamental experiments on the 
interaction of matter with quantized radiation. We refer to the end of this introduction 
for more information and references concerning this aspect of our work. 

According to the fundamental principles of quantum mechanics, states of the sys- 
tems S and £ m are given by normalized vectors (or density matrices) on Hilbert spaces 
TLs and He rn , respectively, [2, 5] 1 . We assume that dimHs < 00, while dimW^ may 
be infinite. Observables As and Ag m of the systems S and £ m are bounded operators 
forming von Neumann algebras Wl$ C BiTLs) an d 2% m C B(H.£ m ). They evolve ac- 
cording to the Heisenberg dynamics R 3 t ^ a l s (As) and R B t 1— > ctg m (A£ m ), where 

X A normalized vector ip defines a "pure" state A t— * {ip,Atp} = Tr(g^A), where = \ip)(ip\. A 
general "mixed" state is given by a density matrix q = J]„>i PnQi> n > where the probabilities p„ > 
sum up to one, and where the ip n are normalized vectors. 
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ag and a £ are *-automorphism groups of 971,5 and 9Jls m , respectively, see e.g. [5]. We 
introduce distinguished reference states, given by vectors ips £ TLs and ip£ m £ 7i-E m - 
Typical choices for ips, ipg m are equilibrium (KMS) states for the dynamics a s , a £m , 
at inverse temperatures (3s, (3g m . The Hilbert space of states of the total system is the 
tensor product 

H = H S ® H c , 

where Tic = ® m >i ^£m> an d where the infinite product is taken with respect to ipe = 
< S> m >i "tPSm- The non-interacting dynamics is the product of the individual dynamics, 
defined on the algebra Wis ®m>i by a s ® m>1 ot E . It is useful to consider the 
dynamics in the Schrddinger picture, i.e. as acting on vectors in Tt. To do this, we first 
implement the dynamics via unitaries, satisfying 

a i # (i # )=e iti #A # e- iti #, tel, andL # ^ # = 0, (1.12) 

for any A# G SPT#, where # stands for either S or £ m . The self-adjoint operators Ls 
and L£ m , called Liouville operators, act on 7is and He m , respectively. The existence 
and uniqueness of L# satisfying (1.12) is well known, under general assumptions on the 
reference states ■0# [5]. In particular, (1.12) holds if the reference states are equilibrium 
states. Let r m > and V m £ Wis <8> 9Jlg m be the interaction time and interaction oper- 
ator associated to S and £ m . We define the (discrete) repeated interaction Schrodinger 
dynamics of a state vector <p G H, for m > 0, by 

U(m)<f> = e" iim • • • e-' lL2 e~' lLl (t>, (1.13) 

where 

Lk = r k L k + r fc ^ L £n 

n^k 

describes the dynamics of the system during the time interval [t± H 1- T k -\,T\ + • • • + 

Tfc), which corresponds to the time-step k of our discrete process. Hence L k is 

L k = L s + L £k + V k , (1.15) 

acting on <S> H.£ k - Of course, we understand that the operator Lg n in (1.14) acts 
nontrivially only on the n-th factor of the Hilbert space Tic of the chain. 

A state g(-) = Tr(p- ) given by density matrix p on H is called a normal state. (A 
density matrix is a self-adjoint, non-negative trace-class operator of unit trace.) Our 
goal is to understand the large-time asymptotics (m — > oo) of expectations 

e (U(m)*OU(m)) = e(aSr(0)), (1.16) 

for normal states £> and certain classes of observables O that we specify below. We 
denote the random repeated interaction dynamics by 

a&(0) = ?7(m)*0?7(m). (1.17) 
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1.3.2 Reduced dynamics and random matrix products 

Let us explain how we link the dynamics to a product of reduced dynamics operators. 
In order not to muddle the essence of our argument, we only consider the expectation 
of an observable As £ QJts, and we take the initial state of the entire system to be 
given by the vector 

ipo = ips®i>c, (1.18) 

where the ips and tpc are the reference states introduced above. (See Section 4 for the 
general case.) The expectation of As at the time-step m is 

(Vo, a^(4s)Vo> = (V>o, Pj~ Ll ■ ■ ■ f Lm A s e~ iZ ™ ■ ■ ■ e'^P^) ■ (1-19) 
We write simply As for A$ ®~Hc, and we have introduced 

m>l 

the orthogonal projection onto TCs'S'C^c- A first important ingredient of our approach 
is to construct operators K k with the properties 

e iZ Me -iZ fe = e iK kAe -iK k ^ ( 121 ) 

K k ip s ®il>C = 0, (1.22) 

where A in (1.21) is any observable of the total system. Equality (1.21) means that 
the operators K k implement the same dynamics as the L k . Relation (1.22) selects 
a unique generator of the dynamics among all that satisfy (1.21). The existence of 
operators K k satisfying (1.21) and (1.22) is linked to the deep Tomita-Takesaki theory 
of von Neumann algebras, c.f. [7] and references therein. It turns out that the K k are 
non-normal operators on TL (while the L k are self-adjoint). We combine (1.21) with 
(1.19) to obtain 

M),aEr(^)Vo> = (^o,Pe iKl ---e iKm PAsi>o) ■ (1-23) 

A second important ingredient of our approach is to realize that 

Pe iKl ■ ■ ■ e iKm P = Pe iKl P ■ ■ ■ Pe lKm P, ( 1 . 24) 

which follows from the independence of the systems £ m (see (4.7)). We identify Pe lKk P 
with an operator M k on Hs, and thus obtain, from (1.23) and (1.24), 

(Y> , a£r(AsM>> = <Y>s, Mi • • • M m A s fo) • (1-25) 

Because the operators My. = Pe lKk P implement a dynamics, we can show (Lemma 4.1) 
that the M& are contractions for some suitable norm ||| • ||| of C d . Moreover, it follows 
from (1.22) that M^ips = ips, f° r an k. The matrices M k satisfy conditions (1) and (2) 
above and are thus RRDOs. 

Remark. 10. The process of reduction of the dynamics explained in this section does 
not involve any approximation. This is in contrast e.g. with master equation techniques, 
where a reduced dynamics is obtained using Born and Markov approximations, or so- 
called weak coupling approximations. Our method allows us to use the structure of 
repeated interaction systems, and to do without all these approximations. 
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1.3.3 Ergodicity of repeated interaction systems 

Let us denote by , uJ G VL W , the random reduced operator process obtained from 
(1.17), (1.25), where the Mj in (1.25) are iid random matrices. We call ot^ the random 
repeated interaction dynamics determined by the RRDO M{uj) = Pe lK ^P. 

Theorem 1.5 (Ergodic theorem for random repeated interaction systems) 

Let ctftf be the random repeated interaction dynamics determined by an RRDO M{u). 
Suppose that p(M(o>) € ■M-(e)) > 0. Then there exists a set C N , s.t. P(f2s) = 1, 
and s.t. for any ZJ £ £1$, any normal state g and any As 6 9Jls, 



where 6 is given by (1.7). 

Remarks. 11. For the normal state g given by the density matrix p = |V ; o)(V'o|> 
(1.26) follows readily from Theorem 1.3 and relation (1.25). To prove the result for 
arbitrary normal states, one uses a property of the reference state ipo called cyclicity 
and separability, [5, 7]. 

12. Our setup allows us to treat systems having various sources of randomness. 
Not only do we allow for random interactions, but also for random characteristics of 
the systems £ m and S (interesting cases are random temperatures and dimensions of 
the £ m and of S). 

13. In [8], we prove a similar result for a general class of so-called instantaneous 
observables (involving also operators on the chain C), and we discuss the physical 
properties of the asymptotic state. See also [7] for the case where M(u>) = M. 



14. As mentioned above, ^Vo> <Xrj {As)"4>oj is fluctuating, it does not converge 

pointwise (in uJ £ ) , but only in the ergodic average sense. 

15. We present in Section 5 the explicit example of a spin-spin repeated interaction 
open quantum system. 

1.4 Related works 

Random matrices. Results on the convergence (of some kind) of products of 
random stochastic matrices, and non-negative matrices, are numerous, see e.g. [14, 
18, 19] and references therein. These results mostly concern properties of the limiting 
distribution, if it exists, in terms of the properties of the distribution of the RRDOs. 
The techniques used to obtain those results rely heavily on the positivity of matrix 
elements. Random matrix products have been heavily studied also for matrices in 
G/rf(IR), see e.g. [9, 15, 1]. Again, the main focus of these works is on the study of 
the properties of the limiting distribution, if it exists, and on the properties of the 
Lyapunov exponents. In this case, the group property of the set of invertible matrices 
G/^(]R) plays a prominent role in the derivation of the results. 

By contrast, besides conditions (1) and (2) defining RRDO, we do not require our 
matrices to be real valued, to have positive entries, or to be invertible. Moreover, we are 




N 



(1.26) 



n=l 
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concerned here with the limiting properties of the products only, not with the limiting 
distribution. On the one hand, we get a.s. convergence results for the products $ n 
(Theorem 1.4). On the other hand, in order to eliminate the unavoidable fluctuations 
in the products \l/ n , we resort to a limit in an average sense, the Cesaro limit (Theorems 
1.2-1.3). Let us also point out that our results show as well that the top Lyapunov 
exponent of products of RRDOs is zero, and that it is almost surely of multiplicity one, 
see Theorem 3.6. 

Getting informations on the fluctuations of the process around its limiting value is 
certainly an interesting and important issue. It amounts to getting informations about 
the law of the vector valued random variable n^ of Theorem 1.4, which is quite difficult 
in general. There are recent partial results about aspects of the law of such random 
vectors in case they are obtained by means of matrices belonging to some subgroups of 
G7d(M) satisfying certain irreducibility conditions, see e.g. [10]. However, these results 
do not apply to our situation. 

Random ergodic theorems for products * n (cj) have been obtained in a more general 
framework in [4]. They prove the almost sure existence of a Cesaro limit for these 
products. However, the identification of the limit is not provided by this general result. 
This identification relies on the detailed properties of the matrices involved, and in 
particular in the separation of a fluctuating part from a decaying part in the dynamical 
process. In this respect, our contribution consists in identifying completely the Cesaro 
limit of products of RRDO. 

Actually, in the case under consideration, Theorem 1.1 yields enough information 
so that Theorem 1.3 becomes a consequence of it and of the existence result proven in 
[4]. Since a direct self-contained proof of Theorem 1.3 is easy to derive, given Theorem 
1.1, we present it for completness. 

We close this comparison by noting that the paper [6] contains results, in a de- 
terministic setup, which are close to some of the results we prove in Section 2. The 
authors of [6] consider infinite products of contractions for some norm on M^(C), in the 
reverse order with respect to RRDP. They prove convergence of the product under the 
following assumptions: i) there exists a subsequence of the matrices appearing in the 
product which converges to a matrix which is paracontracting; ii) the set of invariant 
vectors of this matrix is contained in the set of invariant vectors of all matrices appear- 
ing in the product. This last assumption implies our condition (2) in a deterministic 
setup. However, the products in [6] are in the reverse order with respect to those we 
address, and they address the deterministic framework only. 

Repeated interaction systems. In the experimental setup of a "One-Atom Maser" 
the system S represents one or several modes of the quantized electromagnetic field in 
a cavity, while the elements £ m describe atoms injected into the cavity, one by one, 
interacting with the radiation while passing through the cavity. After interaction, the 
atoms encode certain properties of the field that can be measured after they exit the 
cavity, [16, 21]. In an idealized model for this process, one assumes that all the elements 
£ m represent a copy of the same, fixed quantum system, and that the interaction is given 
by a fixed interaction time and a fixed interaction operator (the same for all m). Such 
idealized repeated interaction systems have been analyzed mathematically in [22, 7]. 
However, it is clear that in actual experiments, neither the interaction time r m (or 
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V m ) nor the elements £ m can be exactly the same for all m\ Typically, the interaction 
time will be random, given e.g. by a Gaussian or by a uniform distribution around a 
mean value, and the state of the incoming atoms will be random as well, for instance 
determined by a temperature that fluctuates slightly around a mean temperature. (In 
experiments, the atoms are ejected from an atom oven, then they are cooled down to a 
wanted temperature before entering the cavity.) It is therefore important to develop a 
theory that allows for random repeated interactions, as we do in the present paper. Our 
approach is based on the techniques of non-equilibrium quantum statistical mechanics 
developed in [7] . We are not aware of any other work dealing with variable and random 
interactions in the quantum context. 

1.5 Organization of this paper 

Our paper is organized as follows. In Section 2, we prove several deterministic results, 
including the fact that * n converges (for fixed w) if and only if there are no fluctuations 
in the matrix product (Proposition 2.5). In Section 3 we prove our main results, 
Theorems 1.1 and 1.2. In Section 4 we explain the link between repeated interaction 
quantum dynamics and matrix products. We apply, in Section 5, our general results 
to an explicit quantum dynamical system, a random repeated interaction spin system. 

2 Deterministic results 

In this section, we derive some algebraic formulae and some uniform bounds which will 
play a crucial role in our analysis. However, there is no probabilistic issue involved 
here and we shall therefore freeze the random variable. We will thus simply denote 
Mj = M(uj), and 

* n :=Mi---M n . (2.1) 
2.1 Decomposition of the Mj 

Let Pi j denote the spectral projection of Mj for the eigenvalue 1 and define 

4>j := Pf^S, Pj:=m^j\. (2.2) 

Note that (ipj\ips) = 1 so that Pj is a projection and, moreover, M^ipj = tpj. We 
introduce the following decomposition of Mj 

Mj := Pj + QjMjQj, with Qj = t- Pj. (2.3) 

We denote the part of Mj in QjC d , by Mq 3 := QjMjQj. It easily follows from these 
definitions that 

PjPk = Pk QjQk = Qj (2-4) 

QjP k = P k Qj =P k - Pj = Qj - Q k . (2.5) 
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Proposition 2.1 For any n, 

*n = \lps)(0n\+M Ql ---M Qn , (2.6) 

where 

6 n = ^ n + M* Q J n _ l + --- + M* Qn ---M* Q ^ 1 (2.7) 

= M*---M 2 Vi (2.8) 

and where {tps^n) = 1- 

Proof. Inserting the decomposition (2.3) into (2.1), and using (2.4), (2.5), we have 

n 

Since Pj = \ips)(ipj\, this proves (2.6) and (2.7). From (2.5), we obtain for any j,k, 

M Qj M Qk = M Q M k = QjMjM k . (2.9) 

Hence, ^ n = Pi Mi ■ ■ ■ M n + Q±M 1 ■ ■ ■ M n = \^s)(M* • • • M^\ + M Ql ■ ■ ■ M Qn , which 
proves (2.8). □ 

The following lemma is useful in the study of the asymptotic behaviour of *$> n . It is 
proven as the passage from (2.7) to (2.8) is, simply by considering the sequence (Mj n ) n 
instead of (M n ) n . 

Lemma 2.2 We have ££ =1 M* ^ • • • M* ^ Vj fc+1 = M* • • • M*J jn+l , for any sequence 
of indices (j n )n- 

2.2 Uniform bounds 

The operators Mj, and hence the product ^ Tt , are contractions on C d for the norm 
1 1 • 1 1 1 . In order to study their asymptotic behaviour, we need some uniform bounds on 
the Pj, Qj, . . . We denote by || • || the Euclidean norm on C d . The operator norm of a 
rank one projection P = \4>){x\ ls simply ||P|| = \\4>\\ ||x|| (see also after (1.1)). Recall 
also that \\ips\\ = 1- 

Lemma 2.3 We have sup {\\Mj n Mj n _ 1 ■ ■ ■ Mj 1 1| , n G N*, j k G N*} = C < oo. In 
particular, ||* n || < Co- 

Proof. This follows from |||Mj||| < 1, j G N*, and from the equivalence of the norms 
|| • || and HI • |||. □ 

As a consequence, we get the following fundamental uniform estimates. 
Proposition 2.4 Let Co be as in Lemma 2.3. Then, the following bounds hold 
1. For any j G N* ; = ||^|| < C and \\Qj\\ < 1 + C . 
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2. sup{||M . n M Q . n _ i ---M Q . i ||, nG N*, j fe eN*}<C (l + C ). 

3. For any n € N* ; ||0 n || < Cg. 

Proof. 1. Let j G N*. By means of Von Neumann's ergodic Theorem, we have 

JV-l 

Pi .- = lim — V Mf. 

1J N->oo N 3 
k=0 

Hence, we get from Lemma 2.3 that \\Pij\\ < Cq. Since \\ips\\ = 1 we also have 
IIV'j'll ^ Co, and hence ||P,-|| = \\if)j\\ < Co and ||Qj|| < 1 + Co- 

2. Using (2.9), we have Mq^Mq^ • • • M Qji = Q jn M Jn M 3n _ x ■ ■ ■ M n , so that 

W M Qi« M Qi n -i ■ ■ ■ M Qn II ^ \\QiA\Co < 0,(1 + Co). 

3. From (2.8) and the above estimates, we obtain \\6 n \\ < Co 1 1 "01 1 1 < Cq- 1=1 

2.3 Asymptotic behaviour 

We now turn to the study of the asymptotic behaviour of ^f n . If M n = M £ Mie) 
(recall the definition of M(e) given before Theorem 1.1), then Mq converges to zero, 
and ^> n converges to the rank one projection P = IV'sKV'I (exponentially fast). The 
following results shows that in the general case, ^> n converges if and only if ip n does. 

Proposition 2.5 Suppose that liim^^ supjUMg^ • • • Mq h ||, jk £ N*} = 0. Then 9 n 
converges if and only if ip n does. If they exist, these two limits coincide, and thus 

lim * n = \ips)(^oo |, 
where ip^ = lim n _ >oo '0 n . Moreover, |V's)(?/'oo| is a projection. 

In general thus, we cannot expect to obtain pointwise convergence of the 6 n , but we 
have to consider the ergodic average of 9 n (uJ) in Theorem 1.2. This is natural in view 
of the interpretation of our results in terms of dynamical systems (a fluctuating system 
does not converge). For instance, if, in the random setting, M(to) = P{oj) = |V'5)(V'( W )I; 
then * n (uJ) = P(uj n ) converges iff tp(u n ) does (this is actually a particular case of 
Proposition 2.5). For a general random vector there is no reason to expect a 

limit, due to the fluctuations of the Nevertheless, Birkhoff's Ergodic Theorem 

asserts that the ergodic limit exists almost surely [1], and that 

1 N 

lim -£P( Wj ) =E(P) = |^>W)| = IV-sXVa a*, 

N^co iv ' — * 

i=i 

where the projection \ips)(ip I is constant in to. Theorem 1.2 gives an analogous result 
in the general case. 
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Proof of Proposition 2.5. Suppose lim0 n = 9oo exists, and let <p n := 6 n — 0qo- We 
have M*^ = M*0 n _! - M* = n - M*<p n _i = ^ + </> n - M* Hence, for 

any n, MqJ)^ = Q*^ + Q;0 n - M* n n _x =: Q*^ + Xn , where, by Proposition 2.4, 
lim^oo X n = 0. For any p G N, we can write (A/$J p 0oo = (MqJP- 1 ^ + (M^J^Xn, 
and thus we obtain, by induction, (Mg n ) p 6' 00 = Q*0oo + Zlfc=o(-^Q„) fc Xn- By assump- 
tion, given any e > 0, we can chose po large enough so that, for any n, ||(-^Q n ) Po || < £■ 
Using Proposition 2.4, we thus get that for any n, 

HQ^ooll < ell^ooll +PoC (Co + l)||Xn||. 

Since x« goes to zero and e is arbitrary, this proves that lim^oo Q^Ooo = 0. Now, 
(ipS, #oo) = 1 because {tps,0 n ) = 1 for all n. Thus we have ip n = P*0oc = &oo ~ Q n ^oo- 
Therefore limV>n exists and equals #oo. 

Conversely, suppose \miip n = ^oo exists. Since M*ip n = ip n for all n, we have 

M*V>oo = M^lpoo - 1p n ) + 1p n . 

Combining this with Lemma 2.3 shows that limM*?/^ = ip^. Furthermore, we have 
for any n,p£N, 

n+p 

e n+p = E M ^w-- M 4 + A- 

3=1 

n n+p 

= M Q n+p --- M Q n+ ,Y, M Q n --- M Qi + ^+ E M L +p --- M Q j+ ^3 

3=1 j=n+l 

= M% n+p ...M% n+1 6 n + M* +p ...M* +1 ip n+1 , 
where we used Lemma 2.2 to obtain the last line. It follows that 

||0n+ p -Y>oo|| < ||M^ +p ---M^ n+1 || \\e n \\ + \\M* +p ---M* +1 ^ n+1 -i; 00 )\\ 

+IIK+ P ---K+iV'oo - V'ooll 

Let e > 0. By assumption, there exists po such that for any n, ||^g n+p ■ ■ ■ ^Q n+1 II < e - 
Thus we have, for all n, 

Pn+po - V'ooll < e||0„|| + C \\ip n +i - V'ooll 

n+po— 1 

+||m: +P0 Voo-^oo+ £ m: +po ---m; +1 (m;^oo-^oo) 

j=n+l 

< eCo + CoIIV'ti+I - V'ooll + \\ M n+p 1poc ~ V'ooll 

+Co(po - 1) sup \\Mjtpoo - V'ooll- 
n+l<j<n+po-l 

Since limM*^/^ = ip^, and since e is arbitrary, we obtain lim# n = ip^. To prove that 
|V's)(V'oo| is a projection, it suffices to note that (V'ruV's) = 1 f° r an n - 1=1 

The previous convergence result relies on the (exponential) decay of the product 
of operators Mq. One can show that such a decay holds, in a deterministic situation, 
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provided there exists 5 > such that ||Mq.|| < e , Vj G N. One may obtain this 
bound for specific models, see for instance Section 5.2. To satisfy this bound, it is not 
only required that the Mj belong to M.ie) f°r an h but also that 

1. we have some uniformity in the spectrum of the Mq . , i.e. there exists a 7 > 0, 
such that 

a(Mj) \ {1} C {z I \z\ < e-T}, Vj G N, 

2. we can prove that the corresponding spectral projectors do not behave too badly 
so that we can still bound Mq j efficiently, i.e. ||(5j||e _7 < e~ s , for all j. 

Such conditions look quite stringent. However, Theorem 1.1 shows that in the random 
setting, a similar exponential decay holds under much weaker assumptions. 

3 Proofs of Theorems 1.1 and 1.2 
3.1 The probabilistic setting 

Let (CI, p) be a probability space underlying the random process. Given a random 
variable f(u), we denote by / := E(/) = / /(cj)dp(a;) its expectation. We define the 
probability measure dP on Q ext := in a standard fashion by 

dP = rij>idpj, where dpj = dp, V? G N*. 

We denote points in Q cxt by Zo. For any r G N* , and for any collection of measurable 
sets in fi, {Aj}i<j< r , we have 

P((u;i,W2, ■ ■ ■ ,U r ) G (Al,A 2 , ■ ■ ■ , A r )) = Tl-LKjKr / dp(uj). 

J A, 

We also define the shift T : Q cxt — > f2 ex t by 

(To;)j = cjj+i, V u; = (wj)j e N G fJ cxt . (3.1) 

T is an ergodic transformation of fi ex t- Further, we assume that the distribution of the 
matrices M is caracterized by a measurable map 

n 3 ui i-> Af(wi) G M d (C). (3.2) 

We write sometimes M{u) instead of M{uj\). Hence, for any set B G M^(C), 

Y>{M(u)eB)= V (M- l (B))= f dp(w), 

JM- 1 {B) 

and similarly for other random variables. With these notations we have for the i.i.d. 

process 

* n (u;): = M(wi)M(w2) ■ ■ ■ M(w n ) (3.3) 
= M(T°L<;)M(r 1 L I j)---M(T n - 1 Lj), 
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where uJ G . 

In the same way as in (2.2), we introduce the random variable G C d defined 

as 

^(wi) := Pi(wi)V 5 , (3.4) 

where Pi(wi) denotes the spectral projection of M(o>i) for the eigenvalue 1, and we 
decompose 

M(wi) := 1^5)^(^1)1 + Mq(wi) = P(wi) + Mq(^) (3.5) 

as in (2.3). Note that ip{oJi) and Mq(u\) define bona fide random variables: u\ 1— > 
Pi(o;i) is measurable since oj\ >— > M(wi) is [3]. Consider finally the process 

0„(w) = M*{T n ' 1 uj)M*{T n - 2 uj)---M*{Tuj)^{uj) (3.6) 

n 
3=1 

Note that 9 n is a Markov process, since 6 n+ i(ZJ) = M* +1 (oj n+ i)9 n (iO) . 
3.2 Decay estimates 

In a first step, we derive a bound similar to the one in (1.5), but only for some "good 
u" . Then we will prove (1.5) via the Borel-Cantelli Lemma. 

Lemma 3.1 Suppose p (M{uj) G ■M(e)) > 0. Then there exists an Mq G M.(e) suc ^ 
that for any e > 0, p (\\M(uj) - M \\ < e) > 0. 

Proof. Let M be such that |||M||| < 1 and Mips = ip S - We decompose M as in (2.3), 
and then M G M( E ) if and only if o-(Mq) C {z G C s.t. | z\ < 1}. Note that it is 
equivalent to sr(Mg) < 1, where sr denotes the spectral radius. For n G N, set 

M {n) := 

|M G M d (C) s.t. |||M||| < l,Mip s = 4>s,o-(M Q ) C {z G C s.t. \z\ < 1 - 1/n}}. 

The sets C M^(C) are compact, and we have M.<e) = Un-^(n)- Thus, since 

p (M(cj) G A1(e)) > 0, there exists no such that p (M(cj) G A4( no )) > 0. The lemma 
then follows with Mq G M( no ) by the following standard argument using the compact- 
ness of M[n ). 

We construct a sequence M/% G M( no ) as follows. For any k G N* we cover .M( no ) 
with a finite number of balls with centres in M( no ) and radius 1/k. For each fixed k, 
since p (M(w) G A^(„ )) > 0, there is at least one of these balls which has a non-zero 
probability. Pick one of these balls and denote its centre by Thus, Vfc G N*, 

p (||M(u;) — Mfc|| < ^) > 0. Since .M( no ) is compact, the sequence Mfc converges to 
some Mq G A4( no ) (up to taking a subsequence). Now, for any given e > 0, take k large 
enough so that 1/k < e/2 and \\M k - M \\ < e/2. Then 

p (\\M(u) - Moll < e) > p (\\M(co) - M k \\ < e/2) > 0. 

□ 
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Lemma 3.2 Suppose p (M(w) & M/e)) > 0- Then, there exists Qi C O suc/i £/ia£ 
p(f2i) > 0, and i/iere exisi constants 5 > and C > suc/t i/iaf Vfc £ N*, 



€ => ||M Q (wi)M Q (u; 2 ) • • • Mq(^)|| < Ce 



-<jfc 



Proof. Let Mq be as in Lemma 3.1 and decompose it as Mq = IV'sKV'ol + Mq . Using 
the Riesz representation, 

for a small contour T encircling the eigenvalue one, one gets that the map M Mq is 
continuous on {M G M^(C) s.t. ||M — Mo|| < n}, for r\ small enough. Thus Lemma 3.1 
shows that p (\\Mq(u) - M Qo \\ < e) > 0. 

Hence, for any e > 0, there exists Q t C SI such that p(fi e ) > 0, and such that for 
any Uj G fi e , we can write 

M Q (a; j ) = Mq + A(^) =: M Qo + A,- 

with ||Aj|| < e. With these notations, we have 

M Q (u 1 )M Q (u; 2 )---M Q {oj k ) (3.7) 

— \ \ M kl ~ 1 A, M k2 ~ kl ~ x \i ■ ■ ■ M ki ~ ki ~ 1 ~ 1 A, M k ~ kj 

j=0 l<fei<fc 2 <---<fcj<A; 

Since Mo G we know that there exist 7 > and C > such that ||Mq o || < 

Ce~ 7n , for any n. Inserting this estimate in (3.7), we get for any k and (uji, ■ ■ ■ ,u>k) G 



n k , 



k 

\M Q ( Ul )M Q ( U2 )---M Q (u k )\\ < E C'+V^-'V 

j=0 l<k 1 <k 2 <—<kj<k 
k 



Ce ~ 7fc S (j) ( Ce7e ) j = Ce" 7fc (l + Ce 7 e)* 



We now choose e small enough so that e~ 7 (l + Ce 7 e) = e~ s for some S > 0, and the 
proof is complete. □ 



3.3 Proof of Theorem 1.1 

Let J7i be as in Lemma 3.2, and consider a product of L matrices 

M Q {u j+1 )M Q (u j+2 ) ■ ■ ■ M Q {u 3+L ). 

The L matrices belong to the set Mq(Oi) with positive probability, and by Lemma 3.2, 
we have \\MQ(uj + \)MQ(uj + 2) ■ ■ ■ MQ{ujj + i)\\ < Ce~ 5L . In other words, for any j G N, 

F(\\M Q (u j+1 )M Q (u j+2 ) ■ ■ ■ M Q (uj j+L )\\ < Ce~ 5L ) =:V> p(^i) L > 0. (3.8) 
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We consider now a sequence of nL matrices and estimate the probability that 

||M Q ( Wl )M Q (u, 2 )---M Q (u; Ln )|| < {Ce~ 5L ) m C m+l . (3.9) 

We decompose the above product into n blocks of length L. Estimate (3.9) holds in 
particular if among these n blocks, there are at least m blocks in which (3.8) holds. 
Indeed, let ji, . . . , j m G {0, ■ ■ ■ , n — 1} be such that 

\\M Q (u; jkL+1 )M Q (u; jkL+2 ) • ■ ■ M Q (u; jkL+L )\\ < Ce~ 5L , k = 1, • • • ,m. (3.10) 

We then isolate the m "good blocks" and write 

\\M Q {u 1 )M Q (u 2 )---M Q {u; Ln )\\ 

m 

< Jl\\M Q (u jkL+1 )---M Q (u jkL+L )\\ \\M Q { m ) ■ ■ ■ M Q {u hL )\\ 

k=l 
m—1 

x II \\M Q (u jkL+L+1 )---M Q (u jk+lL )\\ \\M Q (u jmL+L+1 )---M Q (u nL )\\. 

k=l 

Using Proposition 2.4 and (3.10), we obtain (3.9) with C > Cq(1 + Co). Hence, 

n , n. 

¥(\\M Q {^)---M Q {u Ln )\\ < (Ce- 5L ) m C m+1 ) > ]T ( n )v j (l-V) n - j . (3.11) 

j=m 

We now choose m = [f3n], for some (5 > to be chosen later, and where [•] denotes the 
integer part. For L large enough we have the following bound, uniform in ft, 

for some 7 > 0. For < [3 < V, the following classical tail estimate follows from 
Hoeffding's inequality [11] 

j=\J3n] W j=0 

> 1 - e -2(nP-[/3n]) 2 /« > 1 _ e "<™ (3 13 ) 

for a = 2 (P - (3) 2 . Let us denote 

G n := \u = {ujjj^ s.t. ||Mq(wi)M (w2) • • • M Q (io Ln )\\ < Ce-"" 3 ^. 

Then, (3.11)-(3.13) show that, for all n, we have P(G„) > 1 - e~ an . Let G c n be the 
complement of G n in ri ex t - The first Borel-Cantelli Lemma asserts that 

nG c n ) < 00 => p( n n u m > n G^) = 0. 

n>0 



7?i(l _t?) 
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Since 

P(G£) < e"™, (3.14) 

we get that occurs infinitely often with probability zero. Therefore, there exists 
fii C Slcxt, of measure one, such that for any u G fii, there exists n(uJ) s.t. 

||Mq(wi)M (w2) • • • M Q (o; Ln )|| < Ce" aLr \ if n > n(uJ) (3.15) 

where a = j(3/L. The result of Theorem 1.1 now follows: fix L large enough and write 
n = mL + p (with < p < L). We conclude from Proposition 2.4 and (3.15) that if 
wG fii (of measure one) and if n > no(uJ) := Ln(u;) + L (so that m > n(uJ)), then 

||M (wi) • • • M Q K)|| < CoCe~ aLm < C Ce ap e~ an < C Ce aL e~ an . 

This completes the proof of Theorem 1.1. □ 

Remarks. 1. We can take no (ZU) := min{m > | VA; > m, \\Mq(uji)Mq(lu2) ■ ■ ■ Mq^^W < 
Ce" *}, so that, defining 5n = |u s.t. HMq^iJMq^) ■ ■ ■ M Q (u n )\\ < Ce""™}, we 

have |uJ | no(uJ) = mj = f\k>m9k H fl^-i- Hence, by the estimate (3.14) and our 
definition of a, (3 and a, we get 



E(e an °) < e am P(^_ 1 ) < e CT / L £ e^'W 1 = e CT / L £ e (7/3-2(P-/3) 2 WL < 0O> 

m>0 m>0 m>0 

(3.16) 

for /3 small enough. 

2. By Proposition 2.4, (1.5) holds for all U G S^i and for all n G N* , up to replacing 
C by Czj = Ce an °W. 

As we mentioned in the introduction, if one can chose ip(tu) = tp one may expect 
a convergence without any ergodic averaging. Indeed, the following result is a direct 
consequence of (2.8) and Theorem 1.1 

Proposition 3.3 Let M{uj) be a random reduced dynamics operator. Suppose that 
p(M(w) G Mie)) > and there exists if) G C d such that M(uj)*ip = ijj for all u>. Then 
there exists a set Q\ C N , and a constant a > 0, s.t. P(fii) = 1 and s.t. for any 
uJ G Oi, there exists Cjj so that for any n G N* ; 



M l{ u)...M n{ u)-^ sm 



This result applies in particular to bistochastic matrices with ips = if) = -^(1, ■ ■ ■ ,1). 
Infinite product of doubly stochastic matrices have been considered in e.g. [20], where 
the author looks for sufficient conditions which ensure that such a product converges to 
the rank one projection \ips)(' l Ps\ - However, only deterministic products are considered. 
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3.4 Proof of Theorem 1.2 



The proof of Theorem 1.2 relies on the estimate (1.5), but not directly on the more 
stringent assumption p (M(io) G -M-{E)) > 0- Consider the condition 

(RE) There exist ficS] Fi *,C,a>0 such that P(O) = 1, and such that for any w G O, 
there exists no(uJ) G N so that \\Mq{u\) ■ ■ ■ Mq(w„)|| < Ce~ an if n > no(uJ), and, 
\\M Q (ui) • • • Mg(w n )|| < C, otherwise. 

Theorem 1.1 and Proposition 2.4 shows that (RE) is a consequence of p [M{u) G M(e)) > 
0. Consequently, Theorem 1.2 follows from Theorem 3.5 below. We denote by Pijj 
the spectral projection of M := E(M(w)) for the eigenvalue 1. 

Remark. Note that M(lo)i/js = Y^S f° r an w implies M^s = ip$. In general, it may 
happen that 1 is a degenerate eigenvalue for M even though it is non-degenerate for 
any M{oj). This will however not be the case here as the following lemma shows. 

Lemma 3.4 Suppose Assumption (RE) holds, then M G M(e)- Moreover, we have 
P[j= \^s}(0\, where 9 satisfies 

oo 

9 = (1 - M^-V = £Mq^ = = Pl^s- (3.17) 

k=0 



Proof. We have M = \ips)(ip\ + Mq. Since M(uj)ip s = and (V>(w)|t/>s) = 1 for 
any uj, we get MqiPs = and (ip\ips) = 1- Thus, for any n G N we have (see also 

(2.6)-(2.7)) [ 

M n = \^ s ){^ + M Q *^ + --- + M* Q n Vl+M Q n . (3.18) 

By independence, we have, for any k G N, M Q = J dPM Q (^i) • • • M Q {u) k ) so that 

||MQ fe || < J dP||M Q (o; 1 )---M Q (a; fe )||. (3.19) 

Due to Assumption (RE) and considerations similar to (3.16), we can bound the right 
hand side by 

P{n (cJ) > k}C + F{n (uJ) < k}C e - ak < Ce~ ak , (3.20) 
for some constant C. Hence we have 

\\M^ k \\ < Ce~ ak , Vk G N. (3.21) 

In particular, the spectrum of Mq lies inside the open complex unit disk. 
This estimate and (3.18) prove lim M = \tps){@\ where 

ra-^oo 

oo 
k=0 

Hence M G M.(e) an d P\ m = IV'sK^I- The last identities in (3.17) are immediate since 
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Remarks. 1. In general ip is not an eigenvector of M*, and 8 ^ ip; see Theorem 5.3 
and the remark thereafter. 



2. Without Assumption (RE), instead of (3.21), we can only prove that \\Mq \\ < 
k 

Co(l + Co) and that \\M \\ < Cq, V/c G N, which is not enough to exclude that 1 G 

Theorem 3.5 Suppose Assumption (RE) holds. Then, there exists a set £Iq C S1 n * 
such that P(f2o) = 1 an d> f or all ZJ G tto, 

1 * 

lim — V6UuJ) = 0. 

n=l 

Proof of Theorem 3.5. Recall that T denotes the shift map on f2 ext (see (3.1)). 
Using (2.7), we get 

N N n-1 

= Y.T. M Q^~^--- M Q^ j+1 ^ TjuJ ) 

n=l n=l j=0 

N N-k 

= E E M* Q {T k+ ^uj) ■ ■ ■ M^T^^Pu). (3.22) 
k=i j=o 

Let us introduce random vectors O^, k = 1, 2, . . . , by 

0(*)(oj) = ^(wi,-,^) (3.23) 

= M£(T fc -^)M£(T fc ^) • • • M£(T^(T° W ) 

= M^K)^- 1 )^!,--- , Wfc _i), 

so that, by (3.22), 

AT AT „ Af-fc 



n=l fc=l j=0 

oo N—k .. oo 

E 0(k) (T>uJ)- =: £g{k,N,<D). 



N 

k=l j=0 k=l 



For each fixed k, by ergodicity, there exists a set fl(k) C r2 N * of probability one, such 
that, for all uj G fyfc), the following limit exists 



N-k + 1 



lim g(k,N,&) = lim- -Vfl W T^ 

3=0 

1 M 

l im __!_y^0(fc)(T^j) =E(^( fe )). 
v/^oo M + 1 ^ v ; v ; 



M^oo M + _ 

J=0 
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Therefore, on the set fioo := Dk^n^^k) of probability one, for any k G N, we have by 
independence of the M{uj), 

lim g(k,N,ZJ) =K(6 (k) ) = (M*:) k - 1 lp. (3.24) 
N—kx ^ 

Also, by Assumption (RE) and Proposition 2.4, on f2o = n Doo, a set of probability 
1, we have the bound ||0(*)(:P5;)|| < C G C e an ^ Ti ^ e' 01 ^-^. Hence, for all N large 
enough, and for all 1 < k < N, by ergodicity 

1 N ~ 1 

\\g(k,N,W)\\ < C C— ^ e^^^e"^- 1 ) < 2C , CE(e mto )e- a(fc - 1) , 

3=0 

which is finite by (3.16). For k > N, ||flr(fe, JV,t^)|| = < 2C CE(e an °)e~ a ^- 1 \ Since 
e -a(fe-i) g / 1 (N), we can apply the Lebesgue Dominated Convergence Theorem in 
(3.24) to conclude that limjv_ +00 jj Yln=i On(w) = a.s. (recall Lemma 3.4). □ 

3.5 The Mean Ergodic Theorem and Lyapunov exponents 

We present in this section some results for products "in reverse order" of the form 
<I> n (uJ) := M{u n ) ■ ■ ■ M(u>i), as well as results related to the Multiplicative Ergodic 
Theorem. The following results are standard, see e.g. [1], Theorem 3.4.1 and Remark 
3.4.10 ii). The limits 

A$ (uJ) = lim ($ n (cJ)*$„(aJ)) 1 / 2 « and A*(uJ) = lim (^(uJ)**^)) 1 / 2 " 

n— >oo n— >oo 

exist almost surely, the top Lyapunov exponent of A$(tJ) and of A^(aJ) are the same, 
called 7i(w); 7i(w) is constant a.s., and so is its multiplicity. It is usually difficult to 
prove that the multiplicity of 71 (uJ) is 1. In the present case however, we can get this 
result. 

Theorem 3.6 Suppose p(M(w) G .M(e)) > 0. T/ien i/tere exisi a > 0, a random 
vector 7/00(0;) and £Iq C roi/i P(f2o) = 1 such that for any uJ G Slo anc ^ n G N* 

$n(") - |^)(r/oo(w)|| < C^e~ an , for some Qj. (3.25) 

As a consequence, for any to G 7i(^) *s °/ multiplicity one. 

Proof. We first decompose $ n (aJ) in the same way as * n (oJ). We have, 

$n(uJ) = |V-S>(»m(w)I + MqK) • • • M (wi), (3.26) 

where 

r/ n (cJ) = ^(a;i) + M Q (a;i)*V'(a;2) + --- + M (a;i)*---M Q (a; n _i)*V'K) (3.27) 
= M(wi)* • • ■ M(a; n ) VK) (3.28) 
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and (tps , 7] n (io)} = 1. Using (3.27) and Lemma 2.2, we have for any n,p G N 
ri n+ i +p (uJ) - rj n+1 (u) 

V 

= Mq^Y ■ ■ ■ MqK)*^M q K +1 )* • • • M Q (uj n+j )*^(uj n+j+1 ) 

= Mq(cji)* • • • M Q (uj n )*M(uj n+1 y • • • M(u n+p y^{u n+p+1 ). 

Thus, Lemma 2.3 and Proposition 2.4 imply 

\\r) n+p+ i(u) -r] n+1 (uJ)\\ < Cl\\MQ(u n )M Q (u n -i) ■ ■ ■ Mq(ui)\\, 

uniformly in p > 0. Using Theorem 1.1, there exists J7o with P(fio) = 1 and a > such 
that for any n,p G N, uJ G Slo, there exists s.t. ||% +p+ i (uJ) — r/ n+ i(u7)|| < CgC^e"™". 
Therefore, %(aJ) = r] 00 {u) + Ojj(e~ na ). Together with (3.26) and Theorem 1.1, this 
proves (3.25). 

Now, we can write 

$ n (S7)*d» n (uJ) = |»7oo(57))<»7oo^)| + 057(e- na ) 

where f)(ZJ) = f]oo (&)/ \\ Voo so that \i](jlj)){f)(Lo)\ is an orthogonal projector. Hence, 
by perturbation theory, the non-negative self-adjoint operator $ n (u;)*$ n (uJ) has a sim- 
ple eigenvalue ||?7oo(w)|| 2 + Ouj(e~ na ) and all its other eigenvalues are non- negative and 
bounded from above by Cue~ na , for some Cjj > 0. Thus the 2n'th root of <I>* (uJ)<J> n (aJ) 
has in the limit n — ► oo, a simple eigenvalue 1, and all its other eigenvalues are bounded 
from above by e~ a ^ 2 . □ 



4 From repeated interactions to matrix products 

In this section, we show how to reduce the evolution of states of a R.I. quantum system 
(for observables on the small system S) to the study of a product of non-identical 
matrices given by a RDO process. 

4.1 The repeated interaction model 

We call the time step m the index m > 1 that labels the interaction step during which 
the system S interacts with the element £ m of the chain C = £\ + £2 + £ 3 + • • • , even- 
though the duration r m of this interaction may vary with m, with the understanding 
that inf mG N* r m > 0. 

The Hilbert space 7is is finite-dimensional. We denote by Lg k the free Liouvillean 
of £k, acting on He k . The invariant reference vector tpg k G He k is cyclic and separating 
for 2% fc , it is independent of time and satisfies L£ k ifj£ k = 0. The interaction operator 
between the mth element of the chain and the reference system at time m is a self- 
adjoint operator V m G 9Jts <8> %Rs m ■ The coupled dynamics on 7is ® TLs m is then 
generated by the interaction Liouvillean L m = L$ + Lg m + V m . Moreover, if J m and 
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A m denote the modular data of the pair (SC% <£> %Jle m ,ips ® ^£ m ) ( see e -g- [5]), we 
further require 

A^V m A- 1 / 2 e<m s ® m £m , Vm > 1. 

The evolution operator acting on H = 7^s <8> in the Schrodinger picture at time step 
m is given by the unitary U(m) = e ~ lLm e~ lLm ^ 1 ■ ■ ■ e~ lLl , where L m is given in (1.14). 
Let us finally define the Liouville operator K m at time m, see [7, 12, 17], by 



K m = - 

Its main dynamical features, 

rj-TmLm A iT m L m 



L s + L £m +V m - J m AU 2 V m A^/ 2 J„ 



K m ip s ® ip £>m = 0, 



-iK„ 



for 4 G £Dls®£Dle, t G 



(4.1) 
(4.2) 



(see also (1.21), (1.22)), are proven to hold by using standard relations of the modular 
data A m , J m , see e.g. [7]. Note also the bound ||e ± * Xm || < e T m\\A^ 2 v m A m 1/2 \\ _ 

We proceed along the same lines as in [7], taking into account that in the present 
setup, V m may vary with m. Let Tic ■= ®k>iHe k , where the infinite tensor product 
is understood with respect to the reference vector ipc '■= ®fc>iV , £ fe - Since L m and L £k 
commute for k / m, we can write successively 



,-iLi _ e -inii e -iriEfc>i L ifc 



g tT mJ2k<m L £k e - tT ™Lm g ""mEk>m ^£fe 



so that we obtain 



I7(m) 



exp 



m J — 1 
i=l fc=l 



e~ iTlLl exp 



-JTl 



Li jj+ 



j = l fc>j 

(4.3) 



The unitaries act trivially on Tfe and satisfy U^ipc = V'Ci Vm G N*. The cor- 
responding Heisenberg evolution of an observable O G SDT = 971,5 ® 971c at time m is 
(1.17). Given a normal state £ on we want to understand the limit of large m of the 
expectation value g{a 1 ^ I {0)) = g(U(m)* O U(m)), for observables on the small system 
S, i.e. O = As <g) lc (more general observables, such as the so-called instantaneous 
observables following the terminology of [7], will be considered in [8]). In other words, 
we wish to adress the long time behaviour of the evolution of the state g defined by 
g o a^j on Wl. 

Any normal state is a convex combination of vector states, i.e., its density matrix 
is given by a convex combination of rank-one projections, p = YlnP'A^n) {4>rL\-, where 
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^2 n Pn = 1> an( ^ where the (ft n are normalized vectors in W. It is then clear that 
convergence of the vector states \(p n )(<p n \ implies convergence of p (see [7] for more 
details). Therefore, we will limit our attention to vector states (4>, ■ 4>) from now on. 
Moreover, because tp '■= 4>s®4>C is cyclic for the commutant SDT', it is enough to consider 
<fi of the form 

<j> = B'4> s ® =: B'tP, (4.4) 

for some 

B' = B' s ®% =1 B' n ® n>N l Sn G5Jf, (4.5) 

with B' s G SDt^, i?^ 6 9Jt^ n (up to a vanishing error as N — > oo, see [7]). Hence, we are 
lead to to consider 

<v, GB'r^/to^v) = <^ ^r&o^jiq)^), (4.6) 

for observables O = As (g> 1c- In (4.6) we use that B' commutes with a^j(0). 



4.2 Observables on the small system 



In this paper we restrict our attention to the Heisenberg evolution of observables As 6 
Wis- That is, we consider 



aSKAs ® 1) = (As ® 1) ^ M 



ir\L\ IT2L2 



(A s <g> 1) e" 



-271 



Li ^+ 



where we made use of the fact that U m acts trivially on Us- Considering (4.6), we 
note that due to the properties of the unitary U + (m) and those of the Liouvillean, 
(4.1)-(4.2), we have 



= U + (m)* e 



iT\L\ JIT2L2 



= U + (m)* e iKl e iK2 ---< 



e 

%K. 



ir m L 

(A s <g> t)ip. 



(As <8> 1) e ~ iTmLm e iT ™-i L ™-i 



-ir\L 



Let us introduce Pjv = Us <8> Igj 



P 



P 



where P, 



IV'ffcX^fJ- Due to ( 4 - 5 ), we nave (il>\(B')*B' = (ip\(B')*B'P N . Moreover, introducing 
the unitary operator 



U+(N) = exp 



N-l N 



= C/+(m)Pv, 



which is independent of m, we can write for m > N, 
(^(B'YB'a^As^m 



= (iP,(B')*B'U + (NyP N e lKl e 



(A s <8> 



= (V>, (B')*B'U + (N)*e iKl ■ ■ ■ e iKN P N e iKN+1 ■ ■ ■ e iKm (A s <g> t)ip). 

At this level, introducing the projector P = IL5 ® \ipc)(^c\, satisfying Pip = ip, we 
observe that P N e iKN + 1 ■ ■ ■ e iKm (A s <3> = Pe iKN + 1 ■ ■ ■ e iKm P(A s I)ip. Thus, we 



24 



can invoke Proposition 4.1 of [7] which states that for any q > 1 and any distinct 
integers n±, . . . , n q , 

p e iK «i e iK « 2 . . . e iK n q p — p e iK ni p e iK„ 2 p p e iK nq p (4-7) 

We note that the proof of this statement does not require that the K n 's be identical as 
operators on Hs ® Tig. Therefore, introducing operators Mj acting on Hs by 

Pe iK iP = Mj®\il)c)(?l>c\, or Mj~Pe iK jP, (4.8) 

we arrive at 

m&)*B'aMAs®m (4-9) 
= (V'KBO'B'^+We^ 1 • • • e iKN PM N+1 M N+2 ■ ■ ■ M m (A s ® 1)Y>). 

Hence, given the repeated interactions dynamics, in order to understand the long time 
behaviour of the evolution of any states goa^j on the set of observables on the reference 
system S, it is necessary (and sufficient) to address the asymptotic behaviour of the 
product of matrices Mj ~ Pe iK ^P of the type ^ m = M±M 2 ■ ■ ■ M m on Hs- 

Finally, we recall the main features of these matrices Mj inherited from those of 
e lK i which shows that they indeed satisfy properties (1) and (2) of an RRDO. The next 
lemma is an easy generalization of Proposition 2.1 in [7]. 

Lemma 4.1 ([7]) With the definition (4-8), we have for any j G N*, Mjips = ips- 
Moreover, to any <j> £ Tig corresponds a unique A 6 VJts such that <p = Alps, and 
|||(^||| := ||-A||b(^ s ) defines a norm on Hs- As an operator on Hs endowed with the 
norm \ \\ ■ \ \\, Mj is a contraction for any j G N*. 

4.3 Proof of Theorem 1.5 

For random repeated interactions, Theorem 1.3 shows that there exists a vector 6 6 Hs 
satisfying (6,4>s) = 1; such that, with probability 1, 

1 M 

J im MT,*™M = Ws){0\=n, (4.10) 

m=l 

where ir = \ips)(9\ is a rank one projector in Hs- Therefore, (4.9) yields 

M 

M^oo Ivl ' — * 
m=l 

= (V, (B')*B'U + (N)*e iKl ■ ■ ■ e iKN P(irA s ® !)</>) 

= (4>,(B'yB'U+(Nye lK i---e tK »\4>) (e,A s ^s)n s 

= (ip\(B')*B'il)) (e\A s ^s)H s 

= \\B^\\ 2 (9\A s ^s)h s , 

where H-BVII is arbitrary close to \\il>\\ = 1, for N large enough, see (4.4). 

In other words, (4.10) implies the existence of an asymptotic ergodic state qe on 
SDTs, defined by qe(As) := (9,Asips)'Hs-> ^As £ %Rs-> such that for any normal state q 
on Wl, liniM^oo J2m=i £( a jR>(4s® 1)) = Qe{As), which is the statement of Theorem 
1.5. . - ^ ^ 
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5 Random Repeated Interaction Spin-Spin system 



We consider the application of our general results to the concrete situation where both 
S and £ are two-level "atoms", and where the interaction induces energy exchange 
processes. This is a particular case of the third example in [7]. We present the following 
results without proofs, a complete analysis of this (and more general) examples will be 
given in [8]. 



5.1 Description of the model 

The observable algebra for S and for £ is 2l 5 = 2bj = M 2 (C). Let E S ,E £ > be the 
"excited" energy level of S and of £, respectively. Accordingly, the Hamiltonians are 
given by 



h s 




E s 



and hs 





E e 



The dynamics are given by ot s (A) = e lths Ae %ths and a E (A) = e lth£ Ae ltfl£ . We choose 
the reference state of £ to be the Gibbs state at inverse temperature (3, i.e., 



Tr(e-^M) 
Tr(e-^£) ' 



(5.1) 



and we choose the reference state for S to be the tracial state, Qo,s(A) = ^Tr(A). The 
interaction operator is defined by 



as <8> a* £ + a* s <8> as, 



where 



1 




and a 



# 





1 



are the annihilation and creation operators respectively, of # = S,£ . The Heisenberg 
dynamics of S coupled to one element £ is given by the *-automorphism group t *— > 
e lth \ Ae~ lthx , A 6 21,5 ® Slg, h\ = hs + hg + Xv, where A is a coupling parameter. 

To find a Hilbert space description of the system, one performs the Gelfand- 
Naimark-Segal construction of (3ls,ws) and (^g,ue), see e.g. [5, 7]. In this repre- 
sentation, the Hilbert spaces are given by Hs = Hs = C 2 ® C 2 , and the vectors 
representing los and us are 



1 



V J=l,2 



tpj and ips 



1 



ipi®(fi+ e- f3E£/2 ip2 ® ifr] , (5.2 



respectively, i.e., we have Q^,#{A) = (A <g> # = /% = 0, #s 

Here, we have set 



0. 



1 





and tf2 = 




1 



We now illustrate our results on this model. 



26 



5.2 Variable or random interaction times 

The situation where r varies is particular in the sense that we can chose ip(to) = ip so 
that we don't need to consider ergodic means, see also Remark 8 after Theorem 1.3. 
We consider both the deterministic and the random situations. In the deterministic 
case, we require that the sequence (r n ) n avoids the set TN, where 

T — 2?r 

"~ y/(E s -E e ) 2 + 4\ 2 ' 

For interaction times in TN, there appear resonance phenomena that require a more 
refined analysis of the dynamics. This has already been noticed in [7]. In the random 
setting, we only require that the probability of the random variable t{oj) to avoid 
TN is non-zero. Let us denote by d(r, TN) := inf^ |r — kT\ the distance between the 

interaction time r and the set TN. 



Theorem 5.1 Let (r n ) n be a sequence of interaction times. Assume that there exists 
e > such that d(r n ,TZ) > e for all n. Then there exists a 7 > such that 



Tr(e-P' h s) 



\i>s){e- p ' hs ®1 iAs|+0(e-^ n ), 



where (3' = and where ips is given in (5.2). Consequently, we have for any initial 
state q, and for any A$ £ Wis, 

Q(a n RI {A s %)) = Qp^As) + 0(e-^), 

where Q/3',s i- s the Gibbs state (5.1) of the small system for the inverse temperature (3'. 

Physically, it is reasonable to assume that the interaction time t(uj) takes values in 
an interval of uncertainty, since during an experiment, the interaction time cannot be 
controlled exactly, but rather fluctuates around some mean value. We then have the 
following result. 

Theorem 5.2 Let Q 3 uj — > t(cj) € be a random variable. We assume that 
p (t(uj) ^ TZ) > 0. Then there exist a set C ft w and a constant 7 > such that 
F(Qq) = 1, and, for all ZJ £ there exists Czj s.t. 



*n(w) 



< Ore" 



-771 



Consequently, we have for any initial state g, and for any As £ Wis 
Q{a n £{A s ® %)) = Qfy^As) + O^e" 7 ") a.e. 
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5.3 Random energies of S 

We consider the situation where the energy £ is a random variable. Set 

2 



eoM := 



(e s - E e (u) - y/(E s - E e {u)f + 4A 2 ) +4AV r v / (^£ £ H) 2 +4A2 
(e s - E £ {uj) - ^{E s -E £ {u)) 2 + A\^ 2 + 4A 2 



We show in [8] that this is an eigenvalue of M(u>), and that 1 is a degenerate eigenvalue 
of M{oj) if and only if 

y/(E s -E £ (u)y+4>? 

(in which case e${uj) = 1). In the statement of the next result, we denote the average 
of a random variable f(u>) by /. 

Theorem 5.3 Let Q 3 lo — ► Eg(uj) € B* + be a random variable. We assume that 



(5.3) 



T/ien t/zere exists a set ^7 C f2 N * suc/i i/iai P(^7) = 1, and suc/i that for any u> € Q7, 



1 N 

Jim T7^^n(^) = l^5><PB® 1 

N^oo iV * — ' 



n=l 



where 



PE 



1 - (1 - eor^l - eo)(l - 2Z-; 5 )j l^)^) 

1 + (1 - e^)- 1 (l-e )(l-2Z-/ iS )j \<p 2 ){w\ 



+ 



Consequently, we have for any initial state g, for any As £ Wis, an d for any uj £ 

1 N 

jfto N £ ^{As ® 1c)) = %5 (^>), 



n=l 



where g^ s is the Gibbs state (5.1) of the small system for the inverse temperature 



:= — log ( 2 1 - (1 - et)-\l - e )(l - 2Z~,] S ) 



- 1 



(5.4) 



Remarks. 1. The expression (5.4) for the "asymptotic temperature" is of course 
also valid in the case of Theorem 5.2. However, in that case, (3'(uj) = [3 1 and (3 = (3' . 

2. In general, $ 7^ (3' which reflects the fact that in Lemma 3.4, we usually have 
6 7^ ip, as mentioned in Remark 2 after the proof of the lemma. 

3. Physically, one can imagine that different kinds of atoms are injected in the 
cavity, or that from time to time some "impurity" occurs. It is therefore reasonable to 



28 



consider a finite probability space f2, so that assumption (5.3) is satisfied provided at 

least one atom has excited energy Eg satisfying r ^^ 5 ^ +4A ^ 

4. Of course, one can combine Theorems 5.2 and 5.3 to obtain a result for systems 
where the interaction times, excitation energies (and other parameters, like tempera- 
tures of the "incoming atoms" ) fluctuate. 

Acknowledgements. We thank the ACI Jeunes Chercheurs "Stochastic modeling 
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